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THE UPSTREAM MAGNETIC FIELD OF COLLISIONLESS GRB SHOCKS: CONSTRAINT BY FERMI-LAT 

OBSERVATIONS 
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ABSTRACT 

Long-lived > 100 MeV emission has been a common feature of most Fermi-LAT detected gamma-ray bursts 
(GRBs), e.g., detected up to - lO^s in long GRBs 080916C and 090902B and - lO^s in short GRB 090510. 
This emission is consistent with being produced by synchrotron emission of electrons accelerated to high 
energy by the relativistic coUisionless shock propagating into the weakly magnetized medium. Here we show 
that this high-energy afterglow emission constrains the preshock magnetic field to satisfy lO'^nJ'/^mG < B„ < 

lO^nQ^^mG, where no is the preshock density in unit of 1 cm""*, more stringent than the previous constraint by 
X-ray afterglow observations on day scale. This suggests that the preshock magnetic field is strongly amplified, 
most likely by the streaming of high energy shock accelerated particles. 

Subject headings: acceleration of particles — magnetic fields — shock waves — gamma-rays: bursts 



1. INTRODUCTION 

Diffusive (Fermi) shock acceleration is believed to play an 
important role in gamma-ray burst (GRB) afterglow model, 
where a shock propagating into the medium accelerates elec- 
trons to high energy and then synchrotron/inverse-Compton 
emission arises. Although the phenomenological afterglow 
model works generally well with observational data, the shock 
physics is not understood from first principle. One of the main 
issues is h ow the magnetic field is amplified (see review by 
IWaxmanll2 006). Many groups have tried to addres s this issue 
using numerical plasma simulations (see,e.g., KesheTetal] 
|2009, and references there in), but because the calculations 
are extremely demanding numerically, the picture of field 
growth is still not clear. Thus using observations to constrain 
the shock physics parameters will be helpful in this issue. Af- 
terglow observations were used to constrain the downstream 
magnetic field and the accelerated electron energy distribu- 
tion, at both high (Waxman 1997a; Freedman & Waxman] 
120011) a nd low (Waxman 1997b.: Eichler & Waxman 20051) 
energy. iLi & Wa xman (2006, hereafter LW06) had also used 
X-ray afterglow observations on day scale to give a lower 
limit to the preshock magnetic field amplitude, which sug- 
gests either the shock propagates into a magnetized wind of 
the GRB progenitor or strong preshock magnetic field ampli- 
fication occurs, most likely due to the streaming of the ac- 
celerated high-energy particles. LW06 further suggested that 
observations in higher energy range and/or in later time may 
post more stringent constraint on the preshock magnetic field. 

Recently the Large Area Telescope (LAT) on board the 
Fermi satellite revealed some new features in the GRB high- 
energy emission. The onset of > 100 MeV emission is de- 
layed compared to MeV emission, and lasts much longer, 
as a decaying power-law, after the MeV emission already 
ends. For examples, in the bright GRBs, the extended > 
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lO^s in long GRBs 
short GRB 090510 
decreases to below 



100 MeV emission was detected up to ^ 
0809 16C a nd 090902B and - lO^s in 

dAbdo et al. 2009a. bll3), until the flux 

the LAT sensitivity. Kumar & Barniol DuranI (l2009allb ) pro- 
posed that the whole high-energy burst, including the prompt 
phase, is produced by synchrotron emission in the external 
forwa rd shock. Other authors (Corsi, Guetta & Piro 20()9l 
Gao et al. 2009; Ghirlanda et al. 2010; Ghise lhni et al.ll20ia 
De Pasquale et al. ,2010; Pandev et al.. 2010.) also found that 



the light curve slope, the spectral index and the flux level 
of the extended high-energy emission are consistent with the 
synchrotron afterglow model. Even though it may not be true 
that the whole burst is dominated by forward shock emission 
(see discussion in ®, for the lO^s-scale emission the forward 
shock emission is still most favored over other possible mod - 
els (see a discussion in the introduction of lWang et al.ll2010l) . 

In this paper we show that if the extended > 100 MeV emis- 
sion is produced by the synchrotron emission from afterglow 
shock accelerated electrons, the lO^s-scale > 100 MeV emis- 
sion gives much more stringent constraint on the preshock 
magnetic field, compared to that by X-ray afterglow obser- 
vations. The paper is organized as follows: in §2 we list the 
adopted preassumptions in our analysis; in §3 we derive the 
maximum synchrotron photon energy limited by the radiative 
(IC and synchrotron/jitter) cooling of the accelerated elec- 
trons; the KN correction to the IC cooling rate is presented 
in §4; then §5 gives the constraint to the upstream magnetic 
field; and the discussion on the result is given in §6. Unless 
stated otherwise, we use the common notation Qx = Q/IQ^ 
and c.g.s. unit. 

2. PREASSUMPTION 

We adopt some preassumptions in our following analysis. 

(1) The > 100 MeV emission at late time, > 10^- lO^s, 
is assumed to be generated by the synchrotron emission in 
post-shock (downstream) region by electrons accelerated in 
the afterglow shock. Compared to other models, external- 
shock synchrotron model provides the most natural explana- 
tion for the properties of the extended emission, supporting 
this preassumption. 

(2) The synchrotron photons are assumed to be IC scat- 
tered at most once, i.e., neglecting the second-order IC scat- 
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tering. This is generally available since the electrons in early 
afterglow are very energetic, their further interactions with 
up-scattered synchrotron photons take place in deep Klein- 
Nishina (KN) regime. 

The other minor preassumptions about the afterglow param- 
eters are below. 

(3) The postshock (downstream) electrons' energy equipar- 
tition parameter is 0.1. This value is inferred from after- 
glow modelling, and is consisten t with the clustering of ex- 
plosion energy (Frail et al. 2001) and X-ray afterglow lumi- 
nosity (Freedman & Waxman 2001). 

(4) The postshock magnetic field energy equipartition pa- 
rameter is in the range of 10"^ ^ S 10"^. In the af- 
terglow modellin g, values of 10"^ 5, 5. 10"^ are usu- 
ally obtained (e.g. JPanaitescu & Kumar 1200 lHHarrison et al.l 
|2b01). Moreover, the X-ray afterglow observations indicate 

Bi, > O.lriQ mG, with n the medium density (LW06). An 
lower limit to the downstream magnetic field strength is the 
shock-compressed field, B'^ ^^^^ = 4rB„. Here F is the Lorentz 
factor of postshock plasma and the prime denotes the quanti- 
ties in the frame of the postshock plasma. Thus the lower 
Umit to eb is eB,comp = B'l^^^^/^lirV^nmpC^ = Bl/lnnmpC^. 
With the lower limit to Z?„ from X-ray observations, one has 



comp 



> 10 



(1) 



(We have also discussed the case of much lower e^, which 
does not change much the final constraints; see With 
these assumed values of eg and e^, we actually also have Eg > 
es, which is consistent with the IC component observed in 
some afterglows (e.g., GRB 000926: iHarrison et al.ll2001h . 

(5) The medium density is assumed to be 10"^ ^ n < 
lO^cm"^. The n value is not sure but this adopted range is con- 
sistent with afterglow mod elling (e.g., Panaitescu & Kunw] 
UPOI; Harri son et al.ll2001h and with our knowledge of the in- 
terstellar medium. 

(6) The postshock injected electrons follow a energy dis- 
tribution of due/d^e oc with the power-law index as- 
sumed to be 2 < p < 2.5. These values are usually im- 
phed by the afterglow obse rvations (e.g., Wax mani 1 1 997at 
iFreedman & Waxman Il200ll and references therein), and are 
consistent with the theoretical value /? w 2.2 derived for 
isotropic diffusion of accelerated particles (in the test particle 
limit) in both numerical calc ulations (e.g.,|Ac hterber g et al. I 
I2OOII) and analytical analysis (iKeshet & Waxman 11200511 

3. DIFFUSIVE SHOCK ACCELERATION AND MAXIMUM 
SYNCHROTRON PHOTON ENERGY 

In the diffusive (Fermi) shock acceleration mechanism, 
high-energy charged particles are scattered by upstream and 
downstream magnetic field back and forth, respectively, 
crossing the shock front multiples times, and gradually gain 
energy in each crossing. The particle acceleration depends 
on the upstream and downstream magnetic fields. We refer 
readers to LW06 for detailed discussion and the formula used 
below. 

For a relativistic shock expanding into the medium, with 
postshock plasma of Lorentz factor F 1 and preshock mag- 
netic field B„, the acceleration time of electrons to energy 
-feineC^ is derived to be 



%meC 



(2) 



where the correction factor g accounts for that the upstream 
electrons are not deflected by an angle of 1 /F just with a 
fraction 1 /F of the Larmor time, because the preshock field 
structure may not be uniform (e.g.. lLemoine & Pelletier 2003} 
iLemoine & Revenu 2006) and that the electrons do not return 
the downstream region once they are deflect ed by 1 /F, be- 
cause the shock velocity is decreasing (e.g., iPiran & Nakarl 
I2OIO1) . A conservative value is adopted as g « 10 (LW06). 
Hereafter, primed variables denote parameter values mea- 
sured in the downstream frame, while non-primed ones in the 
upstream frame. 

3.1. IC cooling 

When an electron resides in the upstream region, it may 
suffer from energy loss due to IC scattering the synchrotron 
photons from the downstream region. LW06 derived the IC 
cooling time of an upstream-residing electron with -fe, mea- 
sured in the downstream region, is 



(3) 



Here u'^^^ is the synchrotron radiation energy density measured 
in the downstream frame. Because the electrons are confined 
in a region with thickness much smaller than the size of the 
shock R, the synchrotron radiation energy densities in the up- 
stream and downstream regions that we concern are similar, 
thus Mpj, can be either of them. 

The IC scattering may occur in KN regime, then u'^^ should 
be considered as the "effective" energy density of scattered 
photons. Let us define the photon frequency 



(4) 



so that the IC scattering between electrons with (downstream- 
frame) Lorentz factor 7(. and photons at > uf^ enters KN 
regime. Except for the spectrum of the seed photons are very 
hard (which is not the case here), an electron of 7^ mainly in- 
teracts with photons below vf^. So if KN effect is important, 
the photon energy density in eq.© should be replaced by that 
below z.™(7,), ^ v^). 

An electron that has been accelerated to •jeineC^ must satisfy 
that its upstream IC cooling time is longer than the accelera- 
tion time, t^. > ffl/F, which gives an upper limit to 7^, 



l'e< 



3eB„r 



1/2 



(5) 



For the following derivation, it is convenient to define the 
"effective" downstream Compton parameter Yx for relevant 
electrons of as 



(6) 



with u'g = B'J/^TT being the energy density of postshock mag- 
netic field, and further replace Mpj, in eq.© with Mph(< i^^'^) = 

' Eq. ^6) is an approximati on, h owever, as shown in the Appendix, consid- 
ering the full KN effect (Eq. IB6I ) only gives a correction within 20% in the 
problem. 
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The characteristic frequency of synchrotron pho- 
tons emitted by electrons with 7^ downstream is 
V = Q3T^'^eB'^/2TTmeC, where it has been taken into 
account that the synchrotron radiation peaks at 0.3 times the 
gyration frequency of relevant electrons. Thus the upper limit 
in eq. (|5]l for 7^' implies also a maximum observed energy of 
synchrotron photons. 



hvob < 7.1- 



,1/8 



5/8,3/8, 



-keV, 



(7) 



where Fob is the effective Compton parameter of electrons 
with 7g|3 that emitting synchrotron photons at I'oh- Here E 
is the total (isotropic-equivalent) kinetic energy of the shock, 
n is the medium density, and the downstream field B'^ is cal- 
culated from equipartition as B'^ = {32neBnmpy^^rc, where 
one still needs dynamical evolution of the shock. In late time 
when the shock follows the Blandford-McKee self-similar 
solution, the Lorentz factor F drops with radius R as T = 
(ME/l^immpC^yl^R-^l^ (iBlandford & McKee l[T976h . Tak- 
ing the equal arrival time surface into account, the relation be- 
tween F and observer time f is f = R/A-T^c ( IWaxman1ll997d) . 

We will calculate later in ^ the value of Fob as function 
of a range of eg value, considering the KN correction for IC 
cooling of electrons. 

3.2. Synchrotron/jitter cooling 

During the upstream residence time, the electrons will 
also suffer energy loss due to synchrotron or jitter radia- 
tion when gyrating or being deflected by the upstream field. 
The energy loss rate of an electron of 7,. = F7^ due to syn- 
chrotron and jitter radiation in the upstream field B„ is £ = 
(4/3)i77-c7^(B^/87r), then the cooling time due to this process 
(measured in the upstream frame) is = ^enigC^ /E, i.e., 



tB- 



6TTmeC 
^Tl'eVBl- 



(8) 



The electron successfully accelerated to jenieC^, again, re- 
quires its energy loss timescale larger than the acceleration 
time (eq.||2l), Ib > ta, which gives another limit to 7e, 



7;< 



67re 

(TTgBuT 



1/2 



(9) 



Therefore another limit to the observed energy of synchrotron 
photons is 



£1/8 1/2 3/8 
rtl'ob < 11 ■ — TTT^TeY. 



Bu.-bglt 



3/8 
3 



(10) 



4. KN CORRECTION TO IC COOLING 



We are interested in those electrons emitting synchrotron 
photons in the LAT range. We should take hvob = 100 MeV 
and then ^'^^ and Fob are the Lorentz factor and effective 
Compton parameter of relevant electrons, respectively. Now 
we derive the value of Fob, considering the KN correction. 
For this purpose we should derive the synchrotron spectral 
distribution and then calculate the energy density of the syn- 
chrotron radiation below v^^. As e^, ^ e^, IC cooling may 
dominate synchrotron cooling. If the IC scattering takes place 
in KN regime, the complexity arises since in this case the 



energy distribution of the postshock electrons and the syn- 
chrotron photons are coupled and affect each other 

The injected electron distribution downstream follows a 
power law cx 7J'' with minimum Lorentz factor 7^ = fp (7^') = 
eefpinip / me)r . Here fp is the ratio between 7^ and the aver- 
age Lorentz factor. As p >2 the electron energy is dominated 
by the low-energy end electrons, thus fp < 1. If the electron 
distribution at the low energy end is an abrupt cutoff at 7,„, 
then fp = (p-2)/(p- 1). A smoother turnover at 7^ will lead 
to fp > (p-2)/{p-l). The synchrotron radiation by electrons 
with 7,'j peaks at frequency I'm = 0.3F7,'^eB^/27rm(,c, i.e.. 



„ - 1 ^ N/ inl6Fl/2^1/2 2 



(11) 



We focus mainly on the point when t w lO^s. Electrons with 
7,', mainly interact with synchrotron photons below frequency 



^f^, where 



,KN 



= 5AE-''^e''^e-^ 



f-3,3/2 

iJp h ■ 



(12) 



At high enough electron energy, the radiative (synchrotron 
and IC) loss time is shorter than the adiabatic cooling time, 
i.e, the energy-loss time due to the expansion of the post- 
shock plasma. The adiabatic loss time, t'^^ = 6/?/13cF 
jGruzinov & Waxman I [l999), is longer than the radiative 
cooling time for electrons with Lorentz factors exceeding 7^ = 
JimeC / Aaju'gil + Yc)t'^^. Here Yc is the Compton F-parameter 
for electrons with 7^ = 7^. The characteristic synchrotron fre- 
quency of photons emitted by electrons with 7^ = 7^, Vc = 
0.3TYj-eB'j2TTm,c, is 



ly, = 2.2 X 10''*£5;'"eB"!2no'f3 ''"(1 +F,.)"''Hz. 



1/2-3/2-1-1/2, 



(13) 



Electrons with 7^ = 7^ mainly interact with synchrotron pho- 
tons below a frequency 



v^"" = 5.1 X lO'^Ejf eB,_2n;'"f3 ''^(1 +F,)Hz. (14) 
So we also have 



,1/2,-1/2. 



^ = 2.3 X \Q''E,^ill^n]l^{\+Y,) 



(15) 



For those electrons with 7^ = that emit synchrotron pho- 
tons at hvah = 100 MeV (i^ob = Q3TY^eB'j/2nmec), the corre- 
sponding KN frequency is 



.o?=4.9xlO'^£;^;/_V,-^/^Hz, 



'54 '=B,-2'3 



(16) 



thus we further have 



,KN 



: 3.9 X \Q-%fe-l[lt-l,fftl'\ (17) 



Synchrotron self-absorption may be important at low fre- 
quencies. Consider the extreme case where all the postshock 
electrons emit at the absorption frequency Vg with the electron 
Lorentz factor being 7/, = {2'n:meCVa/Q3TeB'^)^l^. The absorp- 
tion coefficient is then « ATne^B'^/2'-^'JmeCVaf' ■ Through 
di^R/T ss 1 we obtain 



(18) 



Note that this should be taken as the upper limit to Va be- 
cause what we assume is the "extreme case" where the self- 
absorption is strongest, and tB and n values have been plugged 
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with the maximum ones. Compared with the KN frequency 
we have ^ I'a, then reach the conclusion that the syn- 
chrotron self-absorption is negligible. 

There are uncertainties in the afterglow model parameters, 
especially for the poorly constrained eg, which is usually cou- 
pled with n in the afterglow modelling. We should scan all 
the possible parameter space. However, when IC cooling is 
important and KN correction is important to IC cooling, the 
electron distribution (after cooling modified) and synchrotron 
spectrum become quite complicated compared to the cases 
when synchrotron cooling dominates or IC scattering t akes 
place in Th omson Hmit (e.g., iNakar et ai.H2009l: IWang et al.l 
l2OO9ll2OT0i) . It is better to pin down some relations between 
the characteristic frequencies in the spectrum, and cancel the 
other cases in the parameter space. 

As we are going to consider the > 100 MeV emis- 
sion of lO^s scale in GRBs 080916C and 090902B (with 
isotropic equivalent gamma-ray energy of order lO^'^erg) 
and of 10-s scale in GRB 090510 {E^ - lO^^erg), we will 
consider only two cases with (E,t) = (10^'*erg,10^s) and 
(lO^^erg, lO^s). In these two cases, with the presumed pa- 
rameter ranges mentioned in ^ we examine the above calcu- 
lations and find the followings are always satisfied. 



,KN 



KN 



and 



(19) 



(20) 



Depending on the relation between Vm and Vc, there are 
two regimes for the bulk postshock injected electrons: "slow 
cooling" regime with f„, < h'c and "fast cooling" regime with 
i^c < v„,. Consider the critical case when v,n = Vc, then the 
synchrotron emission is peaking at v„,{= Vc). Given > v„, 
the IC scattering does not suffer KN suppression, and if only 

single IC scattering is considered, then Y,„ = = (ee/^fi)'''^ 
jSari & Esin 1120011) . which substituted into eq. dTSl l and ly,,, = 
Vc leads to a critical value 



e«c,-=1.7x 10" 



EiA,e'^_JlnQ 



(21) 



For f.B < ^Bfii- electrons are slow cooling, and vice verse. As 
the wide parameter space may allow both regimes to happen, 
in what follows we derive Fob in these two regimes separately. 

4.1. Slow cooling regime 

In this case the j//^ spectrum of the synchrotron radiation 
peaks at ff, i.e., the total energy density in synchrotron pho- 



tons IS u 



ph.syn 



Mp},(< Vc). Let us first derive 7^. 



As only the electrons with 7^ > 7^ efficiently cool, let us de- 
note 77 = {vm/ i^cf'''^^^'^ < 1 the fraction of postshock injected 
electron energy that is rapidly radiated (Sari & Esin 2001). If 
p m2 then the value of 77 is usually order of unity. A lower 
limit can be obtained by taking Yc = 0, then eqs. ([Til l and ( fT3] l 
implies ri > 0.4(£544,-5eLi/p«o/f3)" ' for p = 2.2. 

As the synchrotron spectrum follows vf^, oc v'^l^ dXv < v^^ 
(X ;y*^"''V2 < < i/^^ and then decreases with 1/ above Vc- 
The spectral shape above Vc is affected b y KN correct i on, bu t 
is not interested to this work. Following IWang et al.l (12010 '), 
we can set up equations for Yc, depending on the relations 



between vf^, v,„ and Vc'. 



Ycil+Yc)-- 



C2 



C3 



Vc < y'c 



< Vc 



(22) 

where, assuming vf^, v ^'/^ for v > Vc, the correction 
factors are, approximately, c\ w |(3-p)(p-2) if v™ < v„, <C 

i^c, C2 ~ p-^ if v,„ ^ vf^ < Vc and C3 « 1 if ^ i^f^- Now 
substituting eqs. ( fTTT l. (T3[ and ( fT4l ) into eq.(l22]l. Yc can be 
solved out. 

Since in the parameter space that we concern, we have con- 
ditions of eqs. ( fT9] l and (|20] i, there are only several interesting 
cases for Fob: 



Yob = Yc < 



V, 



KN 

ob 

KN 



< Vc 



C5 'yr<'^n 



(3-p)/2 



^ob 
,KN 



< Vf^ < Vc, 



C6 Kh < < < K 



where the correction factors are approximately C4 



v^^ < v^'^ < v,„ < Vc, C5 ; 



l(3-p)ifv^^^ <v,„ 
U3 - p)(p - 2) if v^ < v„, « Vc « vf 



KN 

c ' 

(23) 
1 if 

< V^^ < Vc, 



'ob 

and C6 ; 

For the case with the lowest eg value allowed, eg = 10"^, 
and with the other parameters being £54 = f3 = no = fe-i = 1, 
the condition for slow cooling regime is marginally satisfied. 
We derive Fob to be, as shown in Appendix, 



0.99 



C3C6 r]E^!^ rifl 

0.1 ,y\,^i\fii^' 



(24) 



where /? = 2 is used when p appears in indices. We neglect 
the redshift effect so far Considering this effect, i.e., Vob 
Vohi^+z) and f f/(l +z), the r.h.s. of eq. (|24li should be 
multiplied by (1 +z)"^'"'. 

4.2. Fast cooling regime 

When €b > eB.a- is satisfied, we have Vc < v^ (7^ < 7^,). This 
means all the postshock injected electron energy is radiated 
rapidly, 77 = 1 . 

Since in the parameter space we are interested we have 
v^ > v„i, the electrons with 7^ = 7,', do not suffer KN sup- 
pression in IC scattering the synchrotron photons, and the syn- 
chrotron vfi, spectrum is peaking at v„,. Now that 7,,, > 7^, we 
have vf^ > v^^ > v,„, i.e., electrons around 7^ even have less 
KN correction in IC scattering synchrotron photons. In this 
case, the electron distribution at 7^ < 7^, < 7,,, still follows the 
result derived in Thomson limit, dnc/djc « %^ ( Sari & Esini 
12001b . The relevant synchrotron spectrum follows vfi, oc v'^^^ 
atv < Vc, oc v^/^ at Vc <v < v„„ and then turnover at v,n, with 
Mpij jyj, w Mpj,(< v„,). Without KN effect, the Compton param- 
eter Y,n of electr ons with 7,'„ is still the same as the solution in 
Thomson limit (ISari & Esin 11200 ll) . 



1/2 



(25) 



Using the spectral form of synchrotron radiation and depend- 
ing on the relations of v^ with Vc and v^, the Compton pa- 
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TABLE 1 

COMPTON y-PARAMETER OF ELECTRONS EMITTING 100-MEV 
SYNCHROTRON PHOTONS 





£54=f3 = l 
n(cm-^) 
10-2 1 102 


£53 = f2 = 1 

n(cm-^) 
10-2 J 


10-2 
10-5 


0.02c7._i 0.06c8,-i 0.06c8,-i 


0.004c7._i 0.03c8._i 
0.02ci4,_i 0.2c36,_i 



Note. — The other parameters are e,, = 0. 1, fp = 1 and p = 2 if p appears 
in exponents . Here c 1 4 = c 1 C4 , C25 = C2 cs , C36 = C3 cg , and c _ 1 = c v / 1 0- ' . 

rameter of electrons with ^'^^ is given by 
^ 4/3 / ^ 1/2 

Y.^. = Y,A )-(u2^^"'' (26) 



where c? « |(/:>-2)/(/?- 1) if < I'c < t^m, and cg 
(/:'-2)/(/?- 1) if < < I'm- The case of < 



ob 



IS 



neglected provided the condition of eq. ( fT9] l. Thus with helps 
of eqs. (fTTT i. (fTlT l, (fT6] l and jZST l, we now can calculate Fob- 
In the case of high value, eg = 10"^, with the other pa- 
rameters being £54 = f3 = «o = ^e.-i = 1, we get, see Appendix, 



Fob « 0-12 



eg 



3/8 
3 



0.2 f 1/8 5/8 1/2 f ■ 
^54 '^B.-l'^e.-l-l P 



(27) 



Here the redshift effect will add a term (1 +z) on the rh.s.. 

We have also calculated the Fob value in a while parameter 
space of es = 10"^- 10"^ and n = 10"^- lO^cm"-' for both cases 
of (£■,?) = (lO'^'^erg, lO^s) and (lO^^erg, lO's), with the results 
presented in Table [T] 

5. CONSTRAINT ON PRESHOCK MAGNETIC FIELD 

For LAT energy range, we should take hi'ob =100 MeV in 
eq.©, which, then, gives a lower limit to the upstream field. 



5„>14 



1/2 5/8.3/8„ , , 



"54 



V 100 MeV 



(l+z)5''**mG. (28) 



Note redshift effect is added hereinafter, and Fob has been de- 
rived in ^ Plugging the Fob values calculated in the cases 
of (£,?) = (lO^'^erg, lO^s) and (lO^^erg, lO^s) into eq.®, it is 
easy to find that the constraint by the latter case is less strin- 
gent, i.e., the value of lower limit to B„ is smaller. Moreover, 
in the former case, we have 

B,(es = 10-^) > 0.5^^«^/^l +z)V«mG, (29) 



B,(eB = 10-2) > 2 



E, 
8ih 



54 

3/8 



'ob 5/8 



/7I/8 0.12 
^54 



nl'\l + zf^^mG. (30) 



The result of more accurate calculation in eg = 10 case, 
Yobiz = 0) = 1.2 (see Appendix), has been used. We see that 
the relatively relaxed lower limit to B„ appears when taking 
the lowest allowed value, eg = 10"^, but the limit is in the or- 
der of mG, and insensitive to the cb value. 

One may argue that the ^ 10"^ assumption is motivated 
by X-ray afterglow observations on day scale, and may not 



hold at early time, ^ lO^s, that is addressed in this paper. In- 
deed, if relaxing the eg > 10"^ assumption, the lower limit to 
B,i does not change much, and only varies slowly with as 

Furthermore, with /ii'ob = 100 MeV, eq.tflOli gives an upper 
limit to the upstream field. 



B«< 110 



1/8 1/2 3/8 
^^54 ^B.-2"{) 

378 

^1^3 



/Wob 



100 MeV 



(1+z)-^Vg. (31) 



It can be seen that taking (E,t) = (10''"^erg, 10"'s) gives more 
stringent constraint than (£',f) = (10^"' erg, lO^s). As we expect 
eB ^ 10~2, then we have the relatively relaxed constraint 

£.1/8 

B,(eB = 10-2) < 102^|-n^,/^l +z)-5/8mG. (32) 

^1^3 

It is interesting to note that the constraints by IC and syn- 
chrotron/jitter cooling together limit the upstream field in a 
closed range that spans about two orders of magnitude. 

6. DISCUSSION 

We have shown that lO^s-scale, > 100 MeV GRB emis- 
sion recently revealed by Fermi-LAT provides stringent con- 
straints on the upstream magnetic field. A lower limit to the 
magnetic field is obtained by requiring the acceleration time 
of electrons producing > 100 MeV synchrotron photons to be 
shorter than their energy-loss time due to IC scattering the af- 
terglow photons. By scanning the possible afterglow param- 
eter space, the lower limit for the magnetic field is given in 
eq. ( |29l ). Interestingly, an upper limit to the magnetic field is 
also obtained by requiring the acceleration time to be shorter 
than the energy-loss time due to synchrotron/jitter radiation 
upstream. Given a maximum equipartition value of magnetic 
field downstream, the upper limit for the field is given in eq. 
( |32] |. Combining both lower and upper limits, the upstream 
magnetic field is limited in a closed range with two orders of 

magnitude uncertainty, lO^ng'^^mG < B„ < lO^nj'/^mG. 

There should be another lower limit to B„ by requiring that 
the acceleration time of 100-MeV emitting electrons is shorter 
than the dynamical time of the afterglow shock. As LW06, 
we neglect this "age limit" because the constraint is much less 
stringent than the one by IC cooling. Even assuming the ^ 10 
GeV photons from the Fermi-LAT GRBs is produced by ex- 
ternal shock synchrotron r adiation, the constrai nt to B„ by this 
argument, B„ > 0.1 mG ( Piran & Nakaiil2010l) . is still much 
less stringent than the cooling limit here. 

The lower limit to upstream field by Fermi-LAT ob- 
servations is larger than the previous constraint using X- 
ray afterglow observations (LW06) by orders of magni- 
tude. It can be seen that this high amplitude field is 
not likely to be provided by a magnetized wind from the 
GRB progenitor (see discussion in LW06), therefore the 
only reasonable origin of this field is due to magnetic field 
amplification upstream, most likely by the s treaming of 
high energy shock accelerated particles (e.g.. Couch et alj 
2008; Medvedev & Zakutnvava 2009; Pelletier et al. 200 J 
Lemo ine & Pelletieil 120101; iNiemiec et all 12010). However, 
compared to the common few /iG-scale interstellar medium 
field, this means an amplification of field amplitude by at least 
3 orders of magnitude, SB /B > 10^, 01 amplification of field 
energy density by > 6 orders of magnitude, 6ub/ub > 10^. 
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Such high contrast amplification suggests that the ampUfica- 
tion may have nothing to do with the prior interstellar medium 
field. This supports a co-evolution picture for magnetic field 
and accelerated particles in the shock. 

The lower limit to the upstream field may also imply a lower 
limit to the downstream one, i.e., the shock-compression 
field. The downstream field limit is, then, eB.min = CB.comp = 
B^/lwnmpC^. The upstream limit assuming e^ass = 10"^ is 
Bu > O.SnQ^'^mG, which however implies a larger downstream 

limit eB.mm = 4 x 10~^«y'''* > e^ ass, thus there is no self- 
consistence in the case of eg = 10"^. The present lower limit 
to the upstream field is mG scale for a wide range of eb, which 
implies shock compressed downstream field of e^.comp ~ 

10-4, 

therefore to be self-consistent it is required that es > lO"*. 
This is consistent with the result s from afterglow mod elling 
(e.g.. iPanaitescu & Kum"^ 120011; [Harrison et al .112001 1) . Be- 
sides, the upper limit to B„ is close to the maximum value 
by equipartition argument, -S^/Stt « nnipC^. So the upstream 
field might be amplified to be close to equipartition, although 
still dominated by preshock matter in energy. 

An important assumption in our analysis is that the long- 
lived, > lO^s, LAT energy range emission is produced by 
synchrotron radiation from electrons accelerated by afterglow 
shocks. We noticed very recently there is a debate on the ex- 
ternal shock interpreta t ion for the LAT detected GeV em ission 
(iPiran & Nakad 120101: iBarniol Duran & Kumad 120101) . We 
note here that the debate is not relevant to o ur assumption. 
The a rgument against external shock model (iPiran & Nakarl 
l2OT0h concerns the very high-energy, > 10 GeV, photon ar- 
rives a t - lO^s, as observed in GRB 090902B (Abd o et al. I 
12009b"). However we consider the lower-energy, > 100 MeV, 
emission of longer timescale, > lO^'s. As the photon index 
is steeper than 2, the emission energy is dominated by low- 
energy, ^100 MeV, photons. It could be that the dominant 
> 100 MeV component of > lO^s scale and the > 10 GeV 
photon of still earlier arrival time have different origins. On 



the other hand, the properties of the > 100 MeV emission, 
i.e., the spectral index, the light curve slope and the flux level, 
are consistent with the external shock model, supporting our 
assumption. We also note that the > 100 MeV flux, locating 
above Vc, is insensitive to the poorly known model parame- 
ters, i.e, n and cb (e.g., Freedman & Waxman 2001), giving 
us confidence on the external shock synchrotron assumption. 

The LAT-detected > 10 GeV photons are almost impos- 
sible to be produced by synchrotron radiation because the 
synchrotron cooling limits the shock acceleration, therefore 
might need a different origin (Li 2010; Piran & Nakar 2010J. 
Moreover, the bright GRBs show temporal variabili ties in 
> log MeV light cur ves, e.g., GRBs 090 510 (Abdo et al. I 
I2009d) and 090902B (lAbdo et al. I l2009bl) . also in contra- 
diction with external forward shock prediction. Within the 
framework of the standard internal shock model, the late, 
large-radius residual collisions may naturally give an inter- 
pretationJJLi & Waxman 2008; Li 2010). 

iKumar & Barniol Duran (2009a b), using the external 
shock model for the LAT-emission, concluded that the post- 
shock field is shock-compressed only from the circum burst 
medium field. This is because they explained the whole burst, 
including the burst phase, by this model, hence need the exter- 
nal shock MeV emission not to exceed the observed level. Our 
result suggests that this might be impossible, i.e, the preshock 
field is almost certainly amplified somehow in the external 
shock interpretation. In this case, the MeV emission from the 
external shock in the prompt burst phase might be naturally 
suppressed because the MeV-synchrotron emitting electrons 
mainly cool by IC scattering the bulk low-energy photons. 
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APPENDIX 
ANALYTICAL DERIVATION 

We present detailed analytical derivation of Fob values here. Consider slow and fast cooling regimes separately. 



Slow cooling 

We discuss first the dependence of the relations between break frequencies on the parameter space. We will use the broken 
power-law approximations, neglecting the correction factors c,'s. We assume p = 2 for simplicity when p appears in exponents. 
Let v,„ = and using eqs. ( flll i. ( fT3l l. (fl4] i and (|22] |. we get the critical e^, 

_ r 1.1 X 10-"?7-V^2£^no'f3' (« > 1-8 X lO-^-^e-ijfE^lt, cnc'; F, > 1) 

- \ 6.1 X (« < 1.8 X 10-V'e-li//£5"if3 cm-3; < 1) ■ ^"^'^ 

This means 

{ > } ^ { > }min[l.l X 10-"?7-Vp£5l«o%"', 6.1 x 10-'^eti/>o'f3']- (A2) 
and, again, using eqs. ( fTTl i. ( fTSl l. ( fT4l i and ( |22] |. we get another critical e^, 

^ ^ r 1.0 X lO-^r'/'e;^_ff/''E^lnf'tl'' (« > 1.8 x 10-\-h;U ffE^lh cm'^ > 1) 
\ 1.8 X lO-%^''n-''' (n < 1.8 X lO-\-h-^_J-%lt3 cm'^; < 1) ' 



This means 



^, « % { ^ } minlLO X 10-V"'«;,'Y'^'"£;J»o"''f . 1^8 x 10-'£,t"»i'"]. (A4) 
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Now let us calculate Fob value in the case with the lowest value allowed, eg = 10 ^, and with the other parameters being 
£54 = f3 = no = Cc.-i = 1- This case marginally satisfies the condition for slow cooling regime, eg < ^B.a, and, from eq. ( IA4l i. 
satisfies Vc < I'f^ as well. Thus from eq. ( l22l i. we have 7^.(1 +7^) = lO'^eg-ie^Lg??, which gives the solution 

Y,= 10W^hlWl',c,. (A5) 

Substitute eq. (IA5b into ly,,, and i^c (eqs. [TT]and[T3]l, and hence into the definition ry = {v„/i'cf-P~^^l'^, then ry can be solved out to 
be 

r, = (0.57£54eB,-5eLi/>of3 (A6) 

For p = 2.2 and 2.4, 77 = 0.94 and 0.87, respectively, close to unity. 
With eqs.( IA5b and ( l23T l. we get, for eg = 10"^ case. 

Fob = 9.9c3C,<^f .-■4^.;^^^/;^/3f3'/l (A7) 



Fflsf cooling 

Again we discuss the parameter space dependence first. Let = Vc and using eqs. ( fT3] l, ( fTSI l and ( |25] l, we get the critical e^, 

= 1 .6 X \Q-\-^!',E-,ln,"hl" , (A8) 

where 

I J } ^ I < } 1 .6 X 10-'*6;!/^£5i«o'^'f3^'- (A9) 

Now consider the case of high eg value, eg = 10"^, with the other parameters being £54 = 13 = no = e^.-i = 1. These parameter 

values obviously satisfies fast cooling condition, > e^ cr, and that < i^^- From eq. (|26] l we write Fob = Y^iv^ /v^y^^cs, 
i.e.. 

Fob = 0.62c8£5"i/'e^%'e;'_f /7'f3'/^ (AlO) 

NUMERICAL CALCULATION FOR eB = 10-5 CASE 

As the lowest case gives the most relaxed lower limit to the upstream field, we carry more detailed numerical calculation 
in this case, in particular we should calculate the whole synchrotron photon spectrum, considering KN effect, for more accurate 
calculation of the IC cooling of electrons. Consider ee = 0.1, es = 10"^, n = Icm"^, E = lO^^erg, and t = lO^s, and consider the 
downstream shocked-plasma frame. Note for simplicity in this section we neglect primes for the quantities in the downstream 
frame. As discussed previously, ^ eg ci- for taken parameters, thus ^ Vc, i.e., the system could be in either slightly fast or 
slow cooling regime, thus we need to consider the possibility of both regimes. 

Consider first slow cooling regime. We assume Vc < I'f^, as analyzed in ^ which can also be checked later (see below). The 
electron distribution at 7,,, < 7 < 7e still follows the injected form cx The distribution beyond jc is strongly affected by 
radiative cooling. By the electron continuity equation, we have 

cx 7"<^"'V7- (Bl) 
The electron energy loss rate might be affected by KN effect, 

7cx72(l+F)cx72[MB + Mph«i^™)]. (B2) 

So the electron distribution around 7 depends on the photon energy distribution around = nieC^/hj, and the electron distri- 
bution and synchrotron photon spectrum are coupled. However, if the photons that affect the electron distribution at a certain 
electron energy range are not those radiated by the same electrons, then the analysis becomes simple. This is the case here. As 
Vc < i^f^, the electrons just above jc have Mph(< z^*^) ~ Mph.syn, independent of 7, thus the electrons follow the form cx 7"(/'+'). 
As 7 increases to 7 > 7e = m^c^ /hvc, the electrons mainly interact with photons at < Vc (where u^, cx ii~^P~^y^) hence 
Mph(< v^) cx 7"(^"^''/2, and cx 7"(^''"''/2. This scaling holds until 7 > 7m = nieC^ /hv,„, beyond which the electrons mainly 
interact with photons at < v,„ where u^, cx hence Mph(< v™) cx 7"^^"', and cx 7"(/'"'/3) go far we actually assume 
Mph(< v^)/ub > 1, this is true even for electrons with 7,„ (see below), but not true for electron with high enough energy, 7 > 7^, 
where Y{jb) = Mph(< i^f^ = nieC^/h-fB) /ub= 1 ■ i^f^ and hence 75 can be solved since the low energy end of the photon spectrum 
is known, cx i^'^^. At 7 > 75, the electrons lie in deep KN regime and the energy loss is dominated by synchrotron radiation, 
so 7 cx 7^ and rij cx 7"''''^''. Given the electron distribution above, the synchrotron photon energy density per unit frequency, for 
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slow cooling regime, can be given as broken power laws, 



Uu OC < 



jy-(p-l)/2 

j^-30;-l)/4 
j^-Cip-A)l6 



V < Vm 
V,„ <V <Vc 
Vc <V <Vc 
Vc<V < Vm 

v,„ <Vb 

Vb < V. 



(B3) 



Here = v^yain.x) and = Vsynijx) with z^syn(7) = 03j^eBij/2TTmeC. Note I'c and are functions of 7(., while iy,„ and can be 
directly calculated for given parameters, and the discussion of is given later 

Next consider the fast cooling case < Vm). As discussed in ^ since > v,n the electrons at 7c < 7 < 7m follow cx 7"^, 
and then (x 7"(^''"'' just above 7m. At 7m < 7 < 7c, electrons mainly cool by photons at Vc < < t'm (where Ui, cx z/"'/^) 
then Mph(< z^™) cx 7"'/^ and n-y cx 7"(/'+'/2) \i gfjjj larger 7's the electrons cool by photons in u^y cx z^'/^ segment, and then 
rij cx 7"(^'"'/3>. Again at large enough 7 for deep KN regime, synchrotron energy loss is dominant then cx 7"''''^'' (7 > 75). We 
can write for fast cooling regime 

1,-1/2 

j^-(2p-l)/4 
,,-(3p-4)/6 



CX < 



-p/2 



V <Vc 

iy,„<iy < i>,„ 

Vc <V <Vb 

Vb < V. 



(B4) 



Note in this case not only Vc and Vc but also vb is function of 7c since peak at Vc (see below). 

The spectral peak at v^^^^ = mm(vm,Vc) can be calculated directly for given parameters, m^,^^ = nePi/,„„rad, where = A-Tn is the 
postshock electron density, P^^^.^ = V^e^Bj / nieC^ is the maximum spectral synchrotron power per unit frequency by one electron, 
and fad = dR/l'iTc is the adiabatic cooling time of postshock plasma. Now discuss vb- Assuming v^^ < mm{v„,,Vc) (checked 

later), v^^ is given by " u,j^^.Jy /v^^i^y'^dv = ub, and hence 7^ = nicC^/hv^^ and vb = t^.syn(7B) are found, vb is a constant in 
slow cooling regime but a function of 7c in fast cooling regime. 

Thus the only unknown value to determine the whole synchrotron spectrum is 7c, relevant to electrons cool significantly in a 
time of tad, i-e., 7(7c) = 7c/4d, which reads 



lc(l+Yc) = 



birmeC 



1 

Yc= — 

Ub 



Uydv. 



(B5) 



Now we carry numerical calculation to find the root 7c of eq ( IB5l l. with helps of eqs ( IB3l l and ( IB4I ). For given parameters the 
root turns out to be 7c = 1.37(1.60) x lO'* for p = 2.2(2.4). Thus 7c and 7m(= 1-07 x 10"*) are close, and the system is marginally in 
slow cooling regime, consistent with the fact that eb < e^.ci . With this result, it can be checked that the assumptions of v^ 



KN 



and Mph(< v„) > Ub (hence 7^ < 7^) are satisfied. 
With 7c known then we have all the break frequencies known and hence the spectrum u,^, thus the Compton Y parameter is 

given by Y(j) = (1 /ub) Jq u^dv. For the observed 100-MeV photons, 70b = 8.3 x 10 , thus we get Fob = 1-37 (independent of 
p). A more accurate calculation is taking the KN cross section, ctkn, and the relativistic energy transfer in scattering into account. 



7(7) = 



8<7tUb 



dvui, 



' (1-M)Vn[i7 
dfi 



(1-/X)] 



1 + - 



(l-M) 



(B6) 



where /i = cosO and 9 is the angle between the interacting electron and photon. This gives Yo\, = 1.18 for p = 2.2 or 2.4. 
results are similar to that using 7(7) = Mph(< i^'™[7])/mb (eqlHl, indicating that the latter is an excellent approximation. 



The 
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